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Abstract 

We discuss, in the context of precision measurement of Am^ and #13, physics capabilities en- 
abled by the recoilless resonant absorption of monochromatic antineutrino beam enhanced by the 
Mossbauer effect recently proposed by Raghavan. Under the assumption of small relative system- 
atic error of a few tenth of percent level between measurement at different detector locations, we 
give analytical and numerical estimates of the sensitivities to Am^ and sin 2 2#i3 . The accuracies 
of determination of them are enormous; The fractional uncertainty in Am^ achievable by 10 point 
measurement is 0.6% (2.4%) for sin 2 26 13 = 0.05, and the uncertainty of sin 2 26 13 is 0.002 (0.008) 
both at la CL with the optimistic (pessimistic) assumption of systematic error of 0.2% (1%). The 
former opens a new possibility of determining the neutrino mass hierarchy by comparing the mea- 
sured value of Am^ with the one by accelerator experiments, while the latter will help resolving 
the 823 octant degeneracy. 
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I. INTRODUCTION 



Recently, an intriguing possibility was suggested by Raghavan that the resonant 
absorption reaction |3[ 

v e + 3 He + orbital e~ ^ 3 H (1) 

with simultaneous capture of an atomic orbital electron can be dramatically enhanced. The 
key idea is to use monochromatic z/ e beam with the energy 18.6 keV from the inverse reaction 
3 H — > z/ e + 3 He + orbital e~, by which the resonance condition is automatically satisfied. (See 
0,0] for earlier suggestions.) He then suggested an experiment to measure #i 3 by utilizing 
the ultra low-energy monochromatic i7 e beam. Though similar to the reactor # 13 experiments 
001 the typical baseline length is of order 10 m due to much lower energy of the beam by 
a factor of ~150, making it doable in the laboratories. The mechanism, in principle, would 
work with more generic setting in which 3 H and 3 He in are replaced by nuclei A(Z) and 
A(Z-l). 

The author of Ref . then went on to even more aggressive proposal of enhancement by 
a factor of ~ 10 11 by embedding both 3 H and 3 He into solids 0,0] by which the broadening 
of the beam due to nuclear recoil is severely suppressed by a mechanism similar to the 
Mossbauer effect 0. Then, the event rate of the #13 experiment is enhanced by the same 
factor, allowing extremely high counting rate. Thanks to the ultimate energy resolution of 
AE U /E U ~ 2 x 1CT 17 enabled by the recoilless mechanism, he was able to propose a table 
top experiment to measure gravitational red shift of neutrinos, a neutrino analogue of the 
Pound- Rebka experiment for photons [10]. 

In this paper, we examine possible physics potentials of the #13 experiment proposed in 
0,0. The characteristic feature of the experiment, which clearly marks the difference from 
the reactor # 13 measurement, is the use of monochromatic beam apart from the shorter 
baseline by a factor of ~150. Then, the most interesting question is how accurately Am^ 
can be determined. Note that even without recoilless setting, the beam energy width is of 
the order of AE V /E V ~ 10" 5 , and it can be ignored for all practical purposes. It is also 
interesting to explore the accuracy of #13 measurement. In addition to possible extremely 
high statistics, the baseline as short as ~10 m should allow us to utilize the settin g of 
continuously movable detector, the one which once proposed in a reactor #i 3 experiment [llj 
but the one that did not survive in the (semi-) final proposal. The method will greatly help 
to reduce the experimental systematic uncertainties of the measurement. 

We will show in our analysis that the accuracies one can achieve for Am^ and #13 determi- 
nation by the recoilless resonant absorption are enormous. At sin 2 2# 13 = 0.05, for example, 
the fractional uncertainty in Am^ determination is 0.6% (2.4%) and the uncertainty of 
sin 2 29 13 is 0.002 (0.008) both at la CL under an optimistic (pessimistic) assumption of 
systematic error of 0.2% (1%). 

What is the scientific merit of such precision measurement of Am| x and #13? With a 
1% level precision of Am^, the method for determining neutrino mass hierarchy by com- 
paring between the two effective Am 2 measured in reactor and accelerator (or atmospheric) 
disappearance measurement 01 3] would work, opening another door for determining the 



neutrino mass hierarchy. It is also proposed 0] that the #23 octant degeneracy can be resolved 
by combining reactor measurement of # 13 with accelerator disappearance (appearance) mea- 
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surement of sin 2 2023 (^23 sin 2 2613)} (See lid [l9| for earlier qualitative suggestions.) The 
results of the recent quantitative analysis 20], however, indicate that the resolving power of 
the method is limited at small #13 primarily because of the uncertainties in reactor measure- 
ment of 013. Therefore, the highly accurate measurement of Am| x and 6^3 which is enabled 
by using the resonant absorption reaction should help resolving the mass hierarchy and the 
2 3 degeneracies. 

In Sec. El we discuss "conceptual design" of the possible experiments. In Sec. 11111 we 
define the statistical method for our analysis. In Sec. IIVI we present numerical analysis of 
the sensitivities of 0x3 an d Am| x measurement. In Sec. El we complement the numerical 
estimate in Sec. II VI bv giving analytic estimate of the sensitivities. In Sec. IVII we give some 
remarks on implications of our results. In Appendix [21 we give a general formula for the 
inverse of the error matrix. 



II. WHICH KIND OF 9 13 EXPERIMENT? 

We give preliminary discussions on which kind of setting is likely to be the best one 
for experiment to measure Am^ and $13 with use of ultra low-energy monochromatic i7 e 
beam. In this section, we rely on the one-mass scale dominant |21( (or the two-flavor) 
approximation of the neutrino oscillation probability P[y e — > i7 e ), though we use the full 
three-flavor expression in our quantitative analysis performed in Sec. II V I It reads 



AE 



(z7 e _ T7 e ) = 1 - sin 2 20 13 sin 2 — f- ) , (2) 



where the neutrino mass squared difference is defined as Am^ = m 2 — m 2 with neutrino 
masses rrii (i = 1 - 3) 2 and L is a distance from a source to a detector. With E v = 18.6 keV, 
the first oscillation maximum (minimum in P(u e — > u e )) is reached at the baseline distance 

£om = 9.2( Am \ . ) 'm. (3) 

V2.5 x 10" 3 eV 2 / V ; 



While the current value of Am| x which comes from the atmospheric [22] and the accelerator 
[23l | measurement has large uncertainties, it should be possible to narrow down the value 
thanks to the ongoing and the forthcoming disappearance measurement by MINOS 24]| and 



1 Here is a short summary for the parameter degeneracy. It is the phenomenon that there exist multiple 
solutions for mixing parameters, #13, 6*23 and 5, for a given set of measurement of disappearance, v e 
and v e appearance probabilities, and the octant ambiguity of 023 is a mong them. The nature of the 
degeneracy may be characterized as the intrinsic degeneracy of #13 and 8 [141] , which is duplicated by the 
unknow n sig n of Arn^j [l5| and the octant ambiguity of $23 for a given sin 2023 |l6|- For an overview, see 
e.g., HQ. 

2 When we speak about discriminating the neutrino mass hierarchy by comparing the two "large" Am 2 
measured in v e disappearance and disappearance channels, one has to be careful about the definition 
of Am 2 which enter into the survival probabilities [T3 | . While keeping this point in mind, we do not try 
to elaborate the expressions of Am 2 in this paper by just writing it as Amjj in v e disappearance channel 



which may be interpreted as Am 2 ff | e in [T3 | . 
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T2K [25[ experiments. Furthermore, the experiment considered in this paper is powerful 
enough to determine both quantities accurately at the same time, if detector locations are 
appropriately chosen. 

The whole discussion of the #13 experiment must be preceded by the test measurement at 
~10 cm or so to verify the principle, namely to demonstrate that the mechanism of resonant 
enhancement proposed in y], 0] is indeed at work. At the same time, the flux times cross 
section must be measured to check the consistency of the Monte Carlo estimate. Then, one 
can go on to the measurement of #13 and Am^, and possibly other quantities. Because of 
the expected high statistics of the experiment it is natural to think about using spectrum 
informations. In the case of monochromatic beam it amounts to consider measurement at 
several different detector locations. 

Let us estimate the event rate. Although the precise rate is hard to estimate, the numbers 
displayed below will give the readers a feeling on what would be the time scale for the 
experiment. The v e flux from 3 H source with strength S MCi due to bound state beta decay 
is given by 



(-*-) (Jl ■ 1 

VlMCi/ V 10 m 



/,- e = l-4xl0 7 — — cm-V 1 (4) 



where the ratio of bound state beta decay to free space decay is taken to be 4.7 x 10 -3 
based on [2^ . The rate of the resonant absorption reaction can be computed by using cross 
section a Tes and number of target atoms Nt as R = NTfu £ cr res . Without the Mossbauer 
enhancement the cross section is estimated to be a res ~ 10 _42 cm 2 0, 0| based on j^f. Then, 
the rate with target mass Mt without neutrino oscillation is given by 

R = 2-4 x 10- 4 ( S ^ T ) ( 2 day- 1 (5) 
VIMCi-kgy \10mJ y K ' 



An improved estimate in [2| entailed a factor of ~ 10 11 enhancement of the cross section by 
the Mossbauer effect after the source and the target are embedded into solids. Assuming 
the enhancement factor, a res ~ 5 x 10 _32 cm 2 and the rate becomes 

(SM \ f L \ ~ 2 

Therefore, one obtains about 1.2 x 10 6 events per day for 1 MCi source and 100 g 3 He target 
at a baseline distance L = 10 m. If the enhancement factor is not reached the running time 
for collecting the same number of events becomes longer accordingly. 

Thus, once the 3 He (and much easier 3 H) implementation into solid is achieved, the event 
rate is sufficient. The real issue for high sensitivity measurement of Am^ and sin 2 28 13 
is whether the produced 3 H can be counted directly without waiting for decaying back to 
3 He by emitting electron. It is because the long lifetime of 12.33 year 0| of 3 H makes it 
impossible to identify which period the decayed 3 H was produced, resulting in the errors 
of the event rate in each detector location. Possibilities of real-time counting and direct 
counting by extracting 3 H atoms are mentioned in Q. In this paper, we assume that at 
least one of such methods works, and it offers opportunity of direct counting of events. Note 
that the detection efficiency need not to be high because of huge number of events. What is 
important is the time-stable counting rate which allows relative systematic errors between 
measurement at different detector locations small enough. 
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III. STATISTICAL METHOD FOR ANALYSIS 



In this section, we define the statistical procedure for our analysis to estimate the sensitiv- 
ities of Awj! and sin 2 2# 13 to be carried out in the following sections. We aim at illuminating 
general properties of the x 2 under the assumption of the small uncorrelated systematic errors 
compared to the correlated ones. 

A. Definition of \ 2 an{ l characteristic properties of errors 

We consider measurement at n different distances L = (i — 1, 2, ...n) from the source. 
Then, the appropriate form of A% 2 which is suited for analytic study |2j| and is simply 
denoted as \ 2 hereafter, is as follows: 



where N is the number of events computed with the values of parameters given by nature, 
and N exp is the one computed with certain trial set of parameters. a c is the systematic 
error common to measurement at n different distances, the correlated error, whereas u usys ^ 
indicate errors that cannot be attributed to a c , the uncorrelated errors. The example of the 
former and the latter errors are as follows: 

• a c (correlated error): Uncertainties in number of target 3 He atoms, errors in counting 
the number of produced tritium nuclei, errors in calculating resonant absorption cross 
section, errors in estimating the efficiency of counting tritium nuclei 

• &usys,i (uncorrelated error) : Possible time dependences of number of decaying tritium 
nuclei and detection efficiency of events 

Since we consider moving detector setting the list of the thinkable uncorrelated systematic 
errors is quite limited. If the near detector with the identical structure with a movable far 
detector exists the error can, in principle, be vanishingly small. One may think of the errors 
of the order of 0.1%-0.3%. It is because the flux times cross section can be monitored in 
real time by a near detector. In fact, the similar values for uncorrelated systematic error are 
adopted in sensitivity estimate of some of the reactor # 13 experiments such as the Braidwood, 
the Daya Bay, and the Angra projects In near future experiments, somewhat 

larger values are taken, 0.6% in Double-Chooz project j32] and 0.35% in KASKA j33|. 

On the other hand, it may not be so easy to control the correlated systematic error cr c . 
The number of 3 H nuclei may be measured when they are implemented into solid. The 
number of target nuclei times the resonant absorption cross section may be measured in a 
research and development stage with a near detector. Therefore, we suspect that the largest 
error may come from uncertainty in counting rate of the produced 3 H nuclei. Of course, 
reliable estimate of systematic errors a c and a usys requires specification of the site to estimate 
the background caused by n 3 H reaction etc. But, it can be experimentally measured by the 
source on and off procedure, as pointed out in [1]. Lacking definitive numbers for o c at the 
moment, we use a tentative value a c = 10% throughout our analysis. We have checked that 
the results barely change even if we use more conservative number a c = 20%. 
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If the direct counting of 3 H atoms does not work, we may have to expect much larger 
systematic errors, because one has to extract event rate at each detector location only by 
fitting the decay curve. In this case, determination of baseline dependent event rates would 
be more and more difficult for larger number of detector locations. Probably, the better 
strategy without the direct counting would be to place multiple identical detectors (or of 
the same structure) at appropriate baseline distances. Even in this case, it is quite possible 
that the uncorrelated systematic error a usys can be controlled to 1% level, as expected in a 
variety of reactor #13 experiments 



B. Approximate form of \ 2 with hierarchy in errors 

By eliminating a through minimization the \ 2 can De written as 

Ax 2 = x T y _1 f 



(8) 



where x is defined as 



x 



N obs_ N ex P N obs_ N exp jyobs _ jyexp 



N exp 



N exp 



Using the general formula given in Appendix El V 1 is given by 



(9) 
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1 + (£2=i A 

\ UK 



(10) 



where a? l4 = ol 



+ 



1 

JET?. 



„ s;/s _, y ..,,.. By construction, the \ 2 depends upon cr U sys,i an d N? ep only 
through this combination. Therefore, the particular case that will be taken in the next 
section, in fact, includes many cases with different event number but with the same a u . 
Under the approximation a\ { -C a 2 , V~ x simplifies; 



1 



E 



k=l 7? 



(11) 



The remarkable feature of (fTT]) is the "scaling behavior" in which x 2 is independent of 
the correlated error a c , and the sensitivity to sin 2 # 13 and Am^ can be made higher as 
the uncorrelated systematic errors as well as the statistical error become smaller. It may 
be counterintuitive because the leading term of the error matrix V is of order cr 2 . (Se e 
Appendix El) It is due to the singular nature of the leading order matrix, as noted in [34| . 



IV. ESTIMATION OF SENSITIVITIES OF Ara^ AND 13 

We now examine the sensitivities of Am^ and sin 2 2#i3 achievable by the recoilless res- 
onant absorption of monochromatic u e enhanced by the Mossbauer effect. The numerical 
estimate of the sensitivities in this section will be followed by the one by the analytic method 
in Sec. fVl 

The setting of movable detector and the expected high statistics of the experiment make 
it possible to consider the situation that an equal number of events are taken in each detector 



6 



location. Of course, the far a detector from the source, the longer an exposure will take. In 
the following analysis, the number of events are assumed to be 10 6 in each detector location. 
Given the rate in ©, and assuming that the direct counting works, it is obtainable in ~10 
days for 100 g 3 He target even if the detector is located at the second oscillation maximum, 
L = 3L m- On the other hand, the number of events 10 6 is sufficient for our purpose because 
it is unlikely that the uncorrelated systematic errors can be made much smaller than 0.1%. 

We take a "common-sense approach" to determine the locations of the detectors and 
postpone the discussion of the optimization problem. We examine the following four types 
of run, Run I, IIA, IIB, and III, for the measurement. 

• Run I: Measurement at 5 detector positions, L = |Lom, f ^om, ^om, |-^om, and §Z/om 
are considered so that A = Am^L / 4E = ir/10, 37r/10, ir/2, 77r/10, 97r/10 are covered. 

• Run IIA: A setting for precision determination of Am^ by measurement at 10 detector 
positions: L — Li (i — 1, ...10) where L i+ i = L i J r |£om and L\ = |Lom so that the 
range A = to tt is covered. 

• Run IIB: A setting for precision determination of Am^ by measurement at 10 detector 
positions: L = Li (i = 1, ...10) where L i+ % = Li + |Lom and L\ = |Lom so that the 
entire period, A = to 2tt, is covered. 

• Run III: A setting of 20 detector positions: L = L { {i = 1, ...20) where L i+1 = 
Li + |Lom and L\ = |Lom (A = — 2tt). It is to check the scaling behavior of the 
sensitivity with respect to errors. 

In the following two subsections IIV Al and IIVBI we examine the cases of the optimistic 
(cusys = 0.2%) and the pessimistic (cr usys = 1%) systematic errors. We stress here that the 
analyses we will present there contain much more general cases. For example, because of the 
scaling behavior discussed in the previous section, the case with N = 10 6 and cr usys = 0.2% 
is equivalent to = 2 x 10 5 and a usys = 0.0%. Similarly, the case with N = 10 6 and 
Cusys = 1% is equivalent to N = 1.33 x 10 4 and a usys = 0.5%. In the last subsection 
IIVC| we give an estimate of the sensitivities using a tentative setting which may be possible 
without direct counting of 3 H atoms. 

A. Case of optimistic systematic error 

We focus in this subsection on the case of optimistic systematic error, from which one 
may obtain some feeling on the ultimate sensitivities achievable by the present method with 
the four Run options described above. As we mentioned earlier, the correlated systematic 
error a c is taken to be a tentative value of 10% throughout our analysis. The uncorrelated 
systematic error a usys , which is assumed to be equal for all detector locations, is taken to 
be 0.2% in this subsection. 

In Fig.[T]we show in sin 2 26*13 — Am| x plane the expected allowed region by Run I, IIA, IIB, 
and III with number of events 10 6 in each location. Throughout the analysis, the true values 
of Am| x is assumed to be Am^ = 2.5 x 10~ 3 eV 2 . The input values of sin 2 2# 13 is taken as 
0.1 and 0.01 in the left and the right panels in Fig.0 respectively. Throughout the numerical 
analyses in this paper, the other oscillation parameters are taken as: Am 2 x = 7.9 x 10 -5 eV 2 , 
sin 2 6 12 = 0.31, and sin 2 6*23 = 0.5. In each panel, the red-solid, the green-dashed, and the 
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FIG. 1: The expected allowed region by Run I, IIA, IIB, and III with number of events 10 6 in 
each location are depicted. The red-solid, the green-dashed, and the blue-dotted lines are for la 
(68.27%), 2a (95.45%), and 3a (99.73%) CL for 2 DOF, respectively. The input values of the 
mixing parameters are marked by asterisks and they are as follows: Amlx = 2.5 x 1(T 3 eV 2 , 
sin 2 2#i3 = 0.1 and 0.01 in the left and the right panels. 
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TABLE I: The expected fractional uncertainty <5(Am 2 )/Am| 1 (0) in % for the optimistic systematic 
error of a usys = 0.2% reachable by the Run I— III defined in the text. The uncertainties are given 
at la (68.27%) CL for 1 DOF, and the numbers in parentheses are the ones at 3a (99.73%) CL for 
1 DOF. In the left, middle, and right columns, the input value of #13 are taken as sin 2 20i3 = 0.1, 
0.05, and 0.01, respectively. 



&usys — 0.2% 


sin 2 26>i3 = 0.1 


sin 2 20i3 = 0.05 


sin 2 26»i3 = 0.01 


Run type 


5(Am 2 )/Am| 1 (0) (in %) at la (3a) CL 


Run I (5 locations) 


0.84 (2.5) 


1.7 (5.0) 


9-6 (til) 


Run IIA (10 locations) 


0.56 (1.7) 


1.2 (3.5) 


6-0 (l}g) 


Run IIB (10 locations) 


0.28 (0.8) 


0.56 (1.6) 


2.8 


Run III (20 locations) 


0.2 (0.56) 


0.4 (1.2) 


2.0 m 



TABLE II: The expected sensitivity to the sin 2 20i3 for the optimistic systematic error of a usys = 
0.2% reachable by the Run I— III defined in the text. The uncertainties <5(sin 2 26*13) are given at 
la (68.27%) CL for 1 DOF, and the numbers in parentheses are the ones at 3a (99.73%) CL for 
1 DOF. In the left, middle, and right columns, the input value of #13 are taken as sin 2 20i3 = 0.1, 
0.05, and 0.01, respectively. 



&usys = 0.2% 


sin 2 20i3 = 0.1 


sin 2 20i3 = 0.05 


sin 2 20i3 = 0.01 


Run type 


sin 2 20i3 at la (3a) CL 


Run I (5 locations) 


0.1 ±0.0026 (0.0078) 


0.05 ±0.0027 (0.0081) 


0.01 ±0.0028 (0.0085) 


Run IIA (10 locations) 


0.1 ±0.0019 (0.0058) 


0.05 ±0.0020 (0.0061) 


0.01 ±0.0021 (0.0064) 


Run IIB (10 locations) 


0.1 ±0.0017 (0.0050) 


0.05 ±0.0018 (0.0053) 


0.01 ±0.0018 (0.0055) 


Run III (20 locations) 


0.1 ±0.0013 (0.0038) 


0.05 ±0.0014 (0.0041) 


0.01 ±0.0014 (0.0042) 



blue-dotted lines are for la (68.27%), 2a (95.45%), and 3a (99.73%) CL for 2 DOF (degrees 
of freedom), respectively. 

To complement Fig. ^ we give in Table the expected sensitivities to Am 2 x at la and 
3a CL (the latter in parentheses) for 1 DOF for Run I- 1 1 1 . They are obtained by optimizing 
sin 2 26*13 m the fit- For relatively large 6* 13 , sin 2 26*i 3 > 0.05 the expected sensitivities to 
Am| x are enormous. For sin 2 26 , 13 = 0.05 the sensitivities are already less than 1% in Run 
IIA, and is about 0.6% in Run IIB both at la CL. The scaling behavior mentioned at the 
end of the previous section is roughly satisfied, as indicated in Table |U (See Sec. |V]for more 
detailed discussions.) For a small value of 6*i 3 , sin 2 26>i 3 = 0.01 the sensitivity to Am^ are 
much worse, as shown in Table E] They are about 6% in Run IIA, and 3% in Run IIB both 
at la CL. If Run III is carried out it can go down to 2%. 

In TablelTTlthe expected sensitivity to sin 2 26*13 at 1°" and 3a CL (the latter in parentheses) 
for 1 DOF are given. The sensitivities to sin 2 2#i3 can be better characterized by 5(sin 2 26*13), 
not its fraction to sin 2 26*i3, as will be understood in our analytic treatment in Sec. By 
Run I one can already achieve the accuracy of <5(sin 2 26?i 3 ) ~ 0.003, and Run IIA or IIB 
reach to 5(sin 2 26 ) i 3 ) ~ 0.002. The effect of measurement at multiple detector locations on 
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improvement of the sensitivity is relatively minor in the case of sensitivities to sin 2 2613. 
This is in sharp contrast to the case of Am^. 
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FIG. 2: The sensitivity limit of #13 by Run I, IIA, IIB, and III with number of events 10 in 
each location are depicted. The red-solid, the green-dashed, and the blue-dotted lines are for lcr 
(68.27%), 2a (95.45%), and 3a (99.73%) CL for 1 DOF, respectively. 

To show the sensitivity limit on 6*13 achievable by the present method, we present in 
Fig. El the excluded regions in sin 2 2#i3 — Am| x space, assuming the case of no depletion of 
u e flux. The four panels in Fig. 121 correspond to Run I, IIA, IIB, and III. In each panel, 
the red-solid, the green-dashed, and the blue-dotted lines are for lcr (68.27%), 2a (95.45%), 
and 3<r (99.73%) CL for 1 DOF, respectively. The sensitivities indicated in Fig. El is quite 
impressive, which reach to sin 2 2# 13 ~ 0.006 at 2a CL even in Run I, and to sin 2 26 , 13 ~ 0.004 
at the same CL in Run IIB. As expected the improvement by adding more detector locations 
is relatively minor. 



B. Case of pessimistic systematic error 

It might be possible that we end up with the error of ~1% due to, e. g., time dependence 
of the source even though the method of movable detector with direct counting of 3 H works. 
In Fig. El we present the similar allowed region in sin 2 2# 13 — Am^ space obtained by the 
same Run I, IIA, IIB, and III with the same number of events of 10 6 in each location but 
with a pessimistic systematic error of cr usys = 1%. At large 9 13 , sin 2 26*13 = 0.1, we still 
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TABLE III: The expected fractional uncertainty (5(Am 2 )/Am| 1 (0) in % for the pessimistic system- 
atic error of a usys = 1% reachable by the Run I— III defined in the text. The uncertainties are given 
at la (68.27%) CL for 1 DOF, and the numbers in parentheses are the ones at 3a (99.73%) CL for 
1 DOF. In the left, middle, and right columns, the input value of #13 are taken as sin 2 20i3 = 0.1, 
0.05, and 0.01, respectively. The column without number represents that no limit is obtained. 



&usys — 1% 


sin 2 26>i3 = 0.1 


sin 2 26>i3 = 0.05 


sin 2 26»i3 = 0.01 


Run type 


5(Am 2 )/Am| 1 (0) (in %) at la (3a) CL 


Run I (5 locations) 


+4.0 /+12\ 
-3.6 l-10j 


+8.0 f+27\ 
-7.2 \~20J 




Run IIA (10 locations) 


2.4 (t 8 7 ;°) 


5-2 (t\t) 


±32 C — ) 


Run IIB (10 locations) 


1.2 (3.5) 


2.4 (+J 7 j) 


tlt(-) 


Run III (20 locations) 


0.8 (2.4) 


1.6 (5.2) 


tli ( - ) 



TABLE IV: The expected sensitivity to the sin 2 20i3 for the pessimistic systematic error of a usys = 
1% reachable by the Run I— III defined in the text. The uncertainties <5(sin 2 26*13) are given at la 
(68.27%) CL for 1 DOF, and the numbers in parentheses are the ones at 3a (99.73%) CL for 1 
DOF. In the left, middle, and right columns, the input value of 6*13 are taken as sin 2 26*13 = 0.1, 
0.05, and 0.01, respectively. The column without number represents that no limit is obtained. 



@usys — 1% 


sin 2 2(9i 3 = 0.1 


sin 2 26>i3 = 0.05 


sin 2 20i3 = 


0.01 


Run type 


sin 2 26»i3 at la {3a) CL 


Run I (5 locations) 


0.1 ±0.011 (0.034) 


0.05 ±0.012 (0.037) 


0.01 i a013 ( 


+0.037N 


Run IIA (10 locations) 


0.1 ±0.009 (0.026) 


0.05 ±0.009 


0.01 i a009 ( 


+0.028 


Run IIB (10 locations) 


0.1 ±0.007 (0.023) 


0.05 ±0.008 (loioli) 


0.01 ±0.008 


'+0.024-j 


Run III (20 locations) 


0.1 ±0.006 (0.017) 


0.05 ±0.006 (0.018) 


0.01 ±0.006 


'+0.019^ 



have reasonable sensitivities to Am^. For Run IIB and III, for example, the sensitivities 
are about 1-2% level for 2 DOF. At sin 2 2^13 = 0.01, however, the sensitivity to Am 2 x is 
lost except for the one at la CL in Run III. It indicates that the value of 0i3 is close to the 
sensitivity limit, and hence we do not place the figure for it though we did for the case of 
optimistic error, Fig. |21 

For more detailed information on sensitivities with the pessimistic systematic error of 
a usys = 1%, we give in Tables ILTT1 and ITVl the sensitivities at la and 3cr CL (the latter in 
parentheses) for 1 DOF to Am^ and sin 2 20i3, respectively. The column without number 
represents that no limit is obtained, in the similar way as seen in the right panels in Fig. 01 
At relatively large $13, sin 2 2^13 = 0.1 and 0.05, the sensitivities to Am^ remain to be 
good, 1.2% and 2.4% at lcr CL for Run IIB. But, the sensitivity quickly drops and is about 
15% in the same Run for sin 2 20i 3 = 0.01. The sensitivity to sin 2 2^13 is still reasonable, 
<J(sin 2 20i 3 ) ~ 0.008 at la CL for Run IIB even at sin 2 26 13 = 0.01. 

For disappearance measurement of P{y e — > i> e ), sin 2 20 13 = 0.01 is a too small value for 
a pessimistic systematic error of a usys = 1% to retain the sensitivity to Am^. Therefore, 
reduction of the uncorrelated systematic error is the mandatory requirement in this method 
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2.6 



2.4 



2.2 



2.6 



2.4 



2.2 



2.1 



2.6 



2.4 



2.2 



2.6 



2.4 



2.2 



\ \Run I 




0.06 0.08 0.1 0.12 0.14 



Run IIA 



0.06 0.08 0.1 0.12 0.14 



Run I IB 



0.06 0.08 0.1 0.12 0.14 



Run III 



0.06 0.08 0.1 0.12 0.14 

sin 2 26i3 



1.5 



Run I 



/ 



0.02 0.04 0.06 0.0S 



Run IIA 



0.02 \0.04 0.06 O.Of 



Run I IB 



0./02 0.04 0.06 0.0! 



Run III 




68.27% C.L. 
95.45% C.L. 
99.73% C.L. 



0.02 0.04 0.06 0.08 

sin 2 26i3 



FIG. 3: The same as in Fig. ^but with the pessimistic systematic error of a usys l%.The input 
values of sin 2 2#i3 is taken as 0.1 and 0.01 in the left and the right panels. 
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TABLE V: The expected fractional uncertainty 5(Am 2 )/Am^ 1 (0) in %, and <5(sin 2 26*13) for the 
pessimistic systematic error of a usys = 1%, 2%, and 3% reachable by Run with 3 detectors as 
defined in the text. The uncertainties are given at la (68.27%) CL for 1 DOF, and the numbers in 
parentheses are the ones at 3a (99.73%) CL for 1 DOF. In the left, middle, and right columns, the 
input value of #13 are taken as sin 2 2^13 = 0.1, 0.05, and 0.01, respectively. The column without 
number represents that no limit is obtained. 



Run 


sin 2 26»i3 = 0.1 


sin 2 26*13 = 0.05 


sin 2 26*13 


= 0.01 


&usys 


5(Am 2 )/Am 2 1 (0) (in %) at la (3a) CL 


1 % 


KS) (ffil)) 


(til) (") 


"(- 


-) 


2 % 


(-) 


(") 


-(- 


-) 


3 % 


(til) (-) 


"(-) 


-(- 


-) 


&usys 


sin 2 26»i3 at la (3a) CL 


1 % 


0.1 ±0.012 (0.035) 


0.05 ±0.012 (0.037) 


0.01 i a013 


^+0.038^ 


2 % 


0.1 ±0.022 (±g*g?) 


0.05 ±0.024 (± a066 ) 


0.01 i a024 


^+0.069^ 


3 % 


0.1 ±0.033 (l™ 89 ) 


0.05 ±0.035 (i : 094 ) 


0.01 i a035 


^+0.098^ 



for accurate measurement of Am^ at small 6*i 3 . 
C. Case without direct detection of 3 H 

Suppose that the direct detection of 3 H in the target is not possible. Then, we may 
have to take the option of multiple detectors with the same structure, giving up the idea of 
movable detector. In this case, most probably, we have to accept a pessimistic value of the 
uncorrelated systematic error of 1-3%. It will cause two important changes in designing the 
experiment. (1) Number of events that can be accumulated in a reasonable time scale would 
be smaller by a factor of ~10 than the case of direct detection. (2) Number of detectors 
that can be prepared by keeping their identity to suppress the uncorrelated systematic errors 
may be limited. Therefore, 3 detector setting, for example, (in addition to a near detector 
which monitors the flux) would be more practical. 

To understand performance of such reduced setting with larger errors, we have carried out 
the similar \ 2 analysis as done in the previous subsections. We take 3 detector setting with 
tentatively determined baselines L = |Lom ; -^om, and 3L Q m, and assume 10 5 events in each 
detectors. We call the setting as Run 0. The three cases of the uncorrelated systematic errors, 
1%, 2%, and 3%, are examined. In Table IV| presented are the expected fractional uncertainty 
<5(Am 2 )/Am 2 i(0) and sin 2 2^13 for Run 0. With 1% of the uncorrelated systematic errors, 
while a sensitivity comparable to Run I is reached for sin 2 2#i 3 , uncertainty of Am 2 x is larger 
by a factor of ~2 compared to Run I. (Note that the baseline settings are not quite optimized 
in Run I.) For the cases of uncorrelated systematic errors of 2% and 3% the uncertainties 
of sin 2 2#i3 get worse by a factor of ~2 and ~3, respectively. The behavior of sensitivities 
to Am^ are more complicated and no numbers are obtained for uncertainties at 3a level 
for most cases. We note that loss of the sensitivities mainly comes from fewer number of 
detectors with larger systematic errors, but not from an order of magnitude smaller number 
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of events. 

The results obtained above indicate that the direct counting, either real time counting 
or an efficient extraction, of the produced 3 H is mandatory to make the type of experiments 
under discussion useful. 



V. ANALYTIC ESTIMATION OF THE SENSITIVITIES 



We complement our numerical analysis of the sensitivities in the previous section by 
presenting analytical treatment of the uncertainties in the Am^ and #13 determination. In 
particular, we derive analytic formulas for the sensitivities of Am^ and sin 2 29 13 under the 
approximation of small uncorrelated systematic error compared to correlated one, a„ < a c 2 . 
In Sec. IIII[ we have argued, assuming feasible direct counting of 3 H atoms, that the hierarchy 
of errors is very likely to hold. 

We restrict ourselves, in consistent with the numerical analysis done in the previous 
section, to the case that an equal number of events are taken in each baseline, N° = N° bs , 
which may be translated into N^ xp = N exp . We also assume, for simplicity, the case of equal 
uncorrelated systematic error in each detector location, cr 2 ^ = o^ sysi + tfW = cr 2 . Under 

these assumptions, V~ l has a simple form 



k- 1 = 4 



I H„ 



n 



/I- 



I 1-1 



_I 1_ 



\ _I _I _I 

> n n n 



-h \ 



(12) 



I — - J 



where H nxn is an n x n matrix whose elements are all unity, fly = 1 for any i and j. 
It indicates again the independence of the \ 2 011 the correlated error a 2 and the "scaling 
behavior" with respect to the uncorrelated systematic error. Then, the x 2 simplifies: 



X 



na: 



E 



06s 



N 



exp 



h3= 



N obs 

3 

prexp 
3 



(13) 



A. Optimal baselines and sensitivities for two detector locations 

Let us start by examining sensitivities for the case of two detector locations. Because of 
a simple setting with monochromatic z7 e beam we can give explicit expression of x 2 i n terms 
of small deviation of the parameters from the true (nature's) values. For this purpose, we 
note that the number of events is given by 

N ex P(N obs ) = Ua res N T TP ee (6 13 , Am 2 31 , L), (14) 

where f„ e denotes the neutrino flux, a res the absorption cross section, Nt the number of 
target nuclei, T the running time, and P ee is a short-hand notation for P{y e — > u e ). In the 
setting in this section, T is adjusted such that an equal number of events is collected at each 
location of the detector. We recall that N obs denotes the event number computed with the 
true value of the parameters, Am^ = Am 2 (0) and sin 2 26*13 — sin 2 26*13(0), whereas N exp 
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denotes the event number computed with possible small deviations <5(Am 2 ) and <5(sin 2 2#i 3 ) 
from the true values of the parameters. Then, i-th component of x vector in @ is given to 
first order in the deviation as 



N: 



obs 



- 1 



+ 



(Li 



sin 



fAm 2 31 (0)L t 



V 4E 



p(0) 
-Tee 

1 9 „ , , /Am 2 (Q)Li 
-sin 2 2^(0) sin ^ 1^ 



<5(sin 2 2fl 13 ) 



x 



5(Am 2 )Li 
2P 



(15) 



where Plv (Z^) = -P e e[$i3(0), Am 2 1 (0), for which we have used the two-flavor expression 
©• 

For simplicity we restrict our discussion in this section to the analysis with single degree 
of freedom. It is a natural setting for estimating ultimate sensitivities; When we discuss 
sensitivity of Am 2 x we optimize x 2 in terms of sin 2 2#i 3 , and vice versa. Or, one can think 
of the situation that, in determination of Am^, sin 2 2^13 is accurately determined by other 
long-baseline accelerator experiments. Under the approximation a 2 <C a 2 and 



using V' 1 in (IT2l . 



r 



is given for small deviations of #13 and Am^ as follows: 



Xe 



[^(sin 2 2g 13 ) 
2a 2 
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Pee (L\ 
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Am^(0)L2 



,(0), 
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(16) 



^Am 2 



X 



fsin 2 26 



13 



Scr 2 

Am§ 1 (0)Li 
2£ 



<5(Am 2 ) 
A^li(0) 



sin 



Am§ 1 (0)Li 
2£ 



Pie 0) (Za) 



Am|jX0)L2 \ ■ / Am§! (0)1,2 
2B y bm ^ 2E 



Pee (L2) 



-, 2 



(17) 



Now, we can address the problem of optimal baseline and estimate the sensitivities of 



sin 2#i3 and Am^ under the approximations stated above. Since sin 2#i3 < 0.1 |25j] it may 



be a reasonable approximation to set PeV(Li) 
of the section. 



1 in the denominator, as we do in the rest 



1. Optimal setting and sensitivity to sin 2 29\^ 

To maximize (j!6)l one should take L\ as short as possible, and L2 at the oscillation 
maximum, the well known feature in the reactor # 13 experiments. Thus, we take L\ = 
and L2 = Lom which makes the square parenthesis in ()16j) unity. Then, one can obtain the 
sensitivity at Nclo~ CL for 2 detector locations as 



5(sin 2 20 13 ) = 2N CL ^= = 2N C J N . (18) 

For a u = 0.2%, 5(sin 2 2# 13 ) = 2.8 x 10~ 3 (8.5 x 10~ 3 ) at la (3a) CL. For a u = 1%, 
5(sin 2 26 ) 13 ) = 0.014 (0.043) at la (3a) CL, which is not so far from the sensitivities quoted 
in the literatures of the reactor # 13 experiments. 
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2. Optimal setting and sensitivity to Am^ 

The optimal baseline setting is quite different for Am| x . We first recall a property of the 
function xsinx] It has the first maximum xsmx = 1.82 at x — 2.02 (L = 0.64Lqm)j has the 
first minimum a; sin a; = —4.81 at x — 4.91 (L = 1.56Lom) ; and then, the second maximum 
xsinx = 7.92 at x — 7.98 (L = 2.54Lom) ; and so on. For simplicity, we restrict ourselves 
to x < 37r, which means L < 3Lqm so that the running time does not blow up. Then, 
the optimal setting is L\ = 0.64Lom and L 2 = 1.56Lom if we restrict to L < 2Lom, and 
L\ = 1.56Lom and L 2 = 2.54Lom if we allow baseline until L < 3Lom- Despite the factor of 
~4 different baseline lengths, we still assume the equal numbers of events at L = 0.64Lom 
and L = 2.54Lom, which implies ~16 times longer exposure time at the latter distance. 

The x 2 is given approximately by 

.2 ... C ^ ( <K Am 31> 



The coefficient c is 5.5 for L < 2Lqm and 20.3 for L < 3Lqm- Then, we obtain the sensitivity 
at Nqlo~ CL for 2 detector locations as 

6{Am 2 31 ) _ 1 Ou_ 

Am 2 3l (0) sin 2 2^ 13 (0) Ci v^' { ] 

Hence, the sensitivity to Am^ depends very sensitively on sin 2 2#i3. 

With a u = 0.2%, §^r|y = 8.5 x 10" 3 at la CL for sin 2 26 13 = 0.1 if we restrict to 

L < 2Lom- If we allow L < 3Lom, the sensitivity becomes better, ^m^^ = 4.4 x 10~ 3 at 
the same CL. If the systematic error is worse, a u = 1%, the sensitivity at sin 2 2# 13 = 0.1 
becomes to ^g%j = 4.3 x 10" 2 and 5 ^ff^ = 2.2 x 10~ 2 at la CL for L < 2L OM and 
L < 3Lqm cases, respectively. 



B. The problem of n detector locations reduces to 2 location case 

We first show that the problem of optimal setting of n detector locations reduces to the 
case of 2 locations under the assumption of equal number of events in each location. To 
indicate the essential point let us first consider a simplified x 2 °f the form x 2 = Z^ij=i( x i ~ 
Xj) 2 and < Xj < 1, which is the essential part of x 2 f° r sin 2 26'i3, the equation (fTBj) . In 
the case of 2 locations the configuration which maximizes the x 2 { n = 2) is xi = and 
x 2 = 1 and Xmaxi 72 = 2) = 1. It is not difficult to observe that in the case of n-locations the 
configuration which maximizes x 2 { n ) is 

• Even n: n=2M; x=0 appears M times, and x = 1 appears M times, X 2 { n )max — M 2 . 

• Odd n: n=2M+l; x=0 appears M+l times, and x = 1 appears M times, or vice versa, 
X 2 {n) max = M(M + 1). 

Thus, the problem of optimal setting with equal number of events at each detector location 
is reduced to the 2 location case. 
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For x 2 f° r A m 3i in (|17|) the situation is slightly different because the function |xsinx| 
increases without limit as x becomes large. Therefore, mathematically speaking, one can 
obtain better and better accuracies as one goes to longer and longer distances in our setting 
of equal number of events at any detector location. 3 But, since we want to remain to a 
reasonable running time, we have restricted our discussions to baselines limited by L < 2Lqm 
or L < 3Lqm m the 2 location case, the restriction which is kept throughout this section. 
Then, one can show that the same result follows for \ 2 f° r Am 2 ^ W^ - Namely, in the case 
of n locations, the highest sensitivity is achieved at the same baselines L\ and L 2 of the 2 
location case; L\ in [|] ( [|] +1 for odd n) times, and L 2 in [^] times, where [ ] implies 
Gauss' symbol. 

The maximal value of x 2 is, therefore, given by 



X 2 ( n )max = (J^j X 2 ( n = fymax] (even n), 



rr 



X 2 (n)max = [ ) X 2 (n = 2) max \ (odd n), (21) 



It can be translated into the uncertainties at Nql CL as 

1? 



5(sin 2 2MW = \/-5(sin 2 2fl 13 )(n = 2), 

n 



8{&ml x )\ 2 {5{A 



^5jJ (n) = V«l^#J (n = 2) ' (22) 

for even n. For odd n, n in ()22|) must be replaced by ^^1. 5(sm 2 26is)(n = 2) and 
Am" (0) ) ( n = 2) are gi ven respectively by (fTSj) and (j20J). Therefore, the sensitivity gradu- 
ally improves as number of runs becomes larger. 4 

At the end of this subsection, we want to note the followings: The reason why we did 
not take these sets of the optimal distances for (9 13 and Am^ obtained in this subsection 
in the numerical analyses in Sec. II VI is that the sensitivity to Am| x is lost if we tune the 
setting optimal for sin 2 26^3, and vice versa. The reason for this is easy to understand; 
At baselines L\ and L 2 which maximizes xl (xlm 2 )' ^L 2 (Xe) vanishes (approximately 
vanishes), because sinxi = sinx2 = at x\ = and x 2 = vr (sin 2 ^ — sm2 if — \ a t %i — f 
and x 2 — ^r). Nonetheless, we will see, in the following subsections, that the sensitivities 
analytically estimated with optimal baseline distances and the numerically calculated ones 
with baselines taken by "common sense" agree reasonably well with each other. 



3 This explains at least partly the reason why the sensitivities to Am|, and sin 2 2#i3 differ in dependence 
on distance from the source to a detector, as indicated in Fig. 3 in [3fi| . 

4 It is the well known feature in the multiple detector setting in the reactor $13 experiments in which one 
obtains better sensitivity as in (|18|l if the the two identical detectors, the near and the far, are each divided 
into small detectors in the same way, if the uncorrelated systematic error a usys is made to be equal with 
that of the original large detector and if the statistical errors are negligible even for divided detectors. 
This point was emphasized by Yasuda |37| . 
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C. Analytic estimation of the sensitivities; sin 2 2^13 

Let us examine the case of a usys = 0.2% and the number of events N = 10 6 which 
was considered in our numerical analysis in Sec. II VI Then, a u = 0.22%. In the case of 5 
locations, n = 5, 5(sin 2 2fl 13 ) = 2.0 x 10~ 3 (6.1 x 10" 3 ) at la (3a) CL. Similarly for 10 and 20 
locations 5(sin 2 2# 13 ) = 1.4 x 10~ 3 (4.2 x 10~ 3 ) and 0.99 x 10~ 3 (3.0 x 10" 3 ), respectively, at 
la (3a) CL. They compare well with the numbers in Table ITT1 though the latter are obtained 
with not-so-tuned baseline settings. Notice that <5(sin 2 26 ) i 3 ) is independent of 6*13 under the 
present approximation of small deviation from the best fit. 

With a usys = 1% the corresponding sensitivities are <5(sin 2 2#i 3 ) = 9.2 x 10~ 3 (2.8 x 10~ 2 ), 
6.4 x 10~ 3 (1.9 x 10" 2 ), and 4.5 x 10" 3 (1.3 x 10" 2 ) for 5, 10, and 20 locations, respectively, 
at la (3a) CL. They are again roughly consistent with the ones in Table llVl 



D. Analytic estimation of the sensitivities; Am^ 

We examine the cases of restriction L < 2Lom and L < 3Lom- In the first case, the 
maximum of the function xsinx is at x = 2.02 (L = 0.64Lom) where xsinx = 1.82, and 
the minimum at x — 4.91 (L = 1.56Lqm) where xsinx = —4.81. In the case of milder 
restriction L < 3Lqm, the maximum of the function xsinx is at x = 7.98 (L = 2.54Lqm) 
where xsinx = 7.92 , and the minimum at x = 4.91 (L = 1.56Lqm) as above. 

In Table EU we give the fractional uncertainties of Am^, 5(Am 2 )/Am 2 1 (0) in % at la 
CL for the optimistic and the pessimistic systematic errors of a usys = 0.2% and 1% (the 
latter in parenthesis), respectively, obtained by using the equations (j2"Uj) and (J2*2~j) . We do 
not show errors at 3a CL because it is obtained simply by multiplying 3. Over-all, the 
analytically estimated uncertainties are in reasonable agreement with those obtained by the 
numerical analysis in Sec. IIV1 Notice that one has to compare the sensitivities of Run I and 
II A with the case of severer restriction L < 2Lqm, and the ones of Run IIB and III with 
the case of milder restriction L < 3L m, because distances beyond 2L OM are involved in 
the latter runs. The fact that our analytical estimates of the errors are smaller than the 
numerical ones by ~30% or so, apart from approximations involved, is consistent with that 
the latter are based on non-optimal baseline distances. It also implies that the baseline 
setting chosen by the "common sense" used in the numerical analysis in Sec. IIVI is not so 
far from the optimal one, indicating that the sensitivities are rather stable against changes 
of baseline setting. 

Based on the numerical and the analytical estimate of the uncertainties of Am 2 x deter- 
mination, we conclude that sensitivities of less than 1% at 90% CL required for resolution 
of mass hierarchy proposed in 0, k| is in reach if the uncorrelated systematic error of 



&usys = 0.2% is realized, and if #i 3 is relatively large, sin 2# 13 > 0.05 in Run IIB. 

We have confined ourselves to the problem of optimal setting of distances under the 
constraint of equal number of events at each location. To minimize running time for a given 
sensitivity, we have to address the problem of optimal detector locations and exposure times 
for a given total running time. It is left for a future study. 
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TABLE VI: The analytically estimated fractional uncertainties of Am^, 5(Am 2 ) / Am^^O) in 
% are given for the optimistic systematic error of a usys = 0.2% and for the pessimistic one (in 
parentheses) of a usys = 1%. The uncertainties are given at la (68.27%) CL for 1 DOF for the 
cases of 5, 10, and 20 detector locations. The upper three lows are for the case of restricted 
baselines, L < 2Lom, whereas the lower three lows are for the cases of somewhat relaxed baseline 
setting, L < 3Lom- In the left, middle, and right columns, the input value of #13 are taken as 
sin 2 26>i3 = 0.1, 0.05, and 0.01, respectively. 



L < 2Lom 


Pusys — 0.2% {a US y S — 1%) 


number of locations 


<5(Am 2 )/Am^(0) (in %) at l^CL 




sin 2 20i3 = 0.1 


sin 2 26»i3 = 0.05 


sin 2 26»i3 = 0.01 


5 locations 


0.61 (2.8) 


1.2 (5.5) 


6.1 (28) 


10 locations 


0.42 (1.9) 


0.85 (3.8) 


4.2 (19) 


20 locations 


0.30 (1.3) 


0.6 (2.7) 


3.0 (13) 


L < 3Lom 


<5(Am 2 )/Am| 1 (0) (in %) at ltrCL 


5 locations 


0.32 (1.5) 


0.62 (2.9) 


3.2 (15) 


10 locations 


0.22 (0.99) 


0.44 (2.0) 


2.2 (9.9) 


20 locations 


0.15 (0.68) 


0.31 (1.4) 


1.5 (6.8) 



VI. CONCLUDING REMARKS 

In this paper, we have explored the potential of high sensitivity measurement of Am^ and 
#13 which is enabled by using the resonant absorption of monochromatic v e beam enhanced 
by the Mossbauer effect. With baseline distances of ~10 m, the movable detector setting 
is certainly possible. Assuming that the direct detection of produced 3 H atom either by 
real time counting or extraction of 3 H atoms works, we have argued that the uncorrelated 
systematic error can be as small as 0.1%-0.3%, if there equipped a near detector with the 
same structure with the far one. It will allow us to determine AtttJi to the accuracies Am|i 
of ~ 0.3 (sin 2 2#i 3 /0.1)- 1 % at la CL (Run IIB, a usys = 0.2%). The error of sin 2 29 13 is 
also small, sin 2 26*13 — 1-8 x 10~ 3 almost independently of 6?i 3 with the same setting. The 
accuracy of the #13 measurement even in Run I, if the systematic error of 0.2% is reached, is 
comparable with that of the next generation accelerator v e appearance experiments [25, 38]. 
It may exceed the accelerator sensitivity if Run IIB is performed. If the systematic error is 
of 1% level, the sensitivity to #13 is similar to the first-stage reactor #13 experiments even in 
Run IIB. 

What is the scientific merit of the precision measurement of Am 2 x and #13? As we have 
already mentioned in Sec. H] the precision measurement of Am 2 x and #13 will have a great 
impact, at least, to the two of the unknowns in the lepton flavor mixing, the neutrino 
mass hierarchy and resolving the #23 octant degene racy . It would be very interesting to 



racy . 

carry out quantitative analyses of such possibilities 39]. We emphasize that such physics 
capabilities can only be made possible by the direct counting of the produced 3 H atoms 
which makes movable detector feasible. We hope that these exciting possibilities stimulate 
further development of the experimental technology toward the goal. 

What is the additional capabilities of the resonant absorption of monochromatic v e beam? 
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With 18.6 keV of neutrino energy the solar oscillation maximum would be reached at 

/ Am 2 \ _1 

-^soiarOM = 290 ( 8xl0 -5e V 2 ) m. Then, it would be worthwhile to explore the possibility 

of precision measurement of 612 and Am^, as was done for the reactor experiments [36]. 
But, in the present case, neither geo-neutrinos nor i7 e flux from nearby reactors (if any) 
contaminate the measurement. In particular, possible movable or multiple baseline set up 
should allow improvement of accuracy of Am^ determination. 

Detection of CP violating effect due to 5 in the v e disappearance measurement requires 
to go down by a factor of O(10~ 6 ) compared to CP conserving terms [40], Unfortunately, 



it would not be in reach despite great potential sensitivities achievable by the resonant 
absorption reaction. 

Finally, the setup of multiple baseline lengths of ~10 m allows a precision test of the pure 
vacuum oscillation hypothesis by observing a sine curve slightly modified by the solar Am^ 
oscillation. It will constrain various possible sub-leading effects such as de-coherence or new 
neutrino interactions to a great precision. 



APPENDIX A: GENERAL FORMULA FOR X 2 FOR n-UNCORRELATED AND 
^-CORRELATED ERRORS 

We consider a general setting with n-uncorrelated and ^-correlated errors, y 2 is given in 
a form as 

X 2 = x T V~ l x (Al) 
where x = (xi, x 2 , x n ) and Xi = -^W — 1. We also introduce a vector notation for 

i 

parameters o„ (p = 1, ..£) for correlated error as o T = (0:1,0:2, ...ctg) for so called the "pull 
type" x 2 SHE!- Then, x 2 can be written (Hf 

X 2 = min 



mm 

a 



(x - Haf D- 1 (x - Ha) + oFD^a 

(o - A~ l H T D- l x) T A (o - A- l H T D- l x) + x T (LT 1 - D' 1 H A' 1 H T D' 1 ) x 

(A2) 

After minimization by ay, V~ l can be written as 

V- 1 = D- 1 - D^HA^^D- 1 . (A3) 

In the above equations, D u is an n x n matrix D u = diag (cr^ 1; • ■ • , a^ n ), D c is an i x I matrix 
D c = diag {(J 2 1, • ■ ■ , cx^), and H is an n x t matrix whose elements are all unity, H Pti = 1 for 
any p — 1, ...I and % — 1, ...n. Finally, the i x I matrix A is defined by 

A = D' 1 + ^D^H. (A4) 

Note that the statistical errors are incorporated in the uncorrelated errors as a\ { = (J 2 ys . ui + 
jjLp. A simple "path- integral" proof is given by [28] that V = D u + HD C H T . 

Here, we are interested in obtaining the explicit form of U _1 in (jA3|) . We first compute 
A~ l . We note that the matrix A given in (jA4J) can be written explicitly as 



.4 



E T ' ( A5 ) 
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T 



( 1 + ei 1 1 

1 1 + e 2 1 

1 1 1 + e 3 

\ 1 1 1 



where e» in (|A6|) is given by 

^cp \; l= \ ui 

It is easy to show that T _1 is given as 

1 



-i 



1 \ 
1 
1 

l + ej 



(p =!,...£). 



(T ) pp 



1 + Er^p et 



(T -1. 



-l 



€ p e 2 
P 



1 +T* i 



P<? 



Cr 



Then, the second term of V~ x in ()A3|) is given by 

1 



(P ^ ?)■ 



tcrml )' 



(2nd tcrm)/*J 2 _2 

ui uj 

1 



^ ^ -^ip(^4 )ppHpj ~\- ^ ^ ffip^A ^jpgHgj 



/-«/ 



L P 



We thus obtain V 1 as 



a 2 - o 2 a 2 ■ 

UI UI UJ 



<7 



cp 



(A6) 



(A7) 



(A8) 
(A9) 



(A10) 



(All) 
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